Quantum paraelectric behavior of pyrochlore PMN 
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Pb1.83Mgo.29Nb1.71Oe. 39 (PMN) crystallizing in a cubic pyrochlore structure exhibits, as the first 
dielectrics with pyrochlore structure, typical feature of quantum paraelectrics - its permittivity 
continuously increases on cooling and levels off below ~ 30 K without any signature of a structural 
phase transition. Broad-band dielectric spectra do not show any dielectric dispersion in the real 
part of permittivity up to 8.8 GHz. THz and infrared spectra reveal a soft polar optic mode which is 
responsible for the temperature dependence of the permittivity. The leveling-off of the permittivity 
at low temperatures obeys the Barrett formula and the fitted vibrational zero-point energy |fcsTi 
corresponds to the measured soft mode frequency. The number of observed infrared phonons exceeds 
that predicted from the factor-group analysis which indicates that the structure is at least locally 
non-cubic. 



PACS numbers: 63.20.-e; 77.22.-d; 78.30.-j; 67.20.+k 
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I. INTRODUCTION 

PbMg 1 / 3 Nb 2 /303 crystallizing in the cubic perovskite 
structure is one of the most studied dielectrics as the 
typical representative of relaxor ferroelectrics. It ex- 
hibits a high (~10 4 ) and strongly diffused peak in the 
temperature dependent permittivity whose position re- 
markably shifts from 245 K (at 100 Hz) to 320 K (at 
1 GHz) Its crystal structure remains cubic down to liq- 
uid He temperatures. The peculiar dielectric properties 
are caused by a wide dielectric relaxation which broadens 
and slows down on coolingji^ The relaxation originates 
from the dynamics of polar clusters, which develop below 
the Burns temperature T^ = 620 The broad distri- 
bution of relaxation frequencies has its origin in random 
fields and random forces as a consequence of chemical 
disorder in the perovskite B sites occupied by Mg 2+ and 
Nb 5+ . The local ferroelectric instability in the polar clus- 
ters was indicated by an unstable polar optic phonon, 
which softens to T<j and hardens above T^il Although 
the dielectric properties of perovskite PbMgi/ 3 Nb 2 /303 
were intensively studied during the last fifty years, not all 
is completely understood, particularly the complex and 
broad dielectric dispersion. 

The related compound Pb1.83Mgo.29Nb!. 71 06. 39 
(PMN) with a pyrochlore structure, which frequently 
grows in the perovskite PMN ceramics and thin films 
as a second phase and significantly deteriorates its 
dielectric properties, has been much less studied. The 
only report known to the authors is that by Shrout and 
Swartz£ They investigated the dielectric response up to 
400 kHz down to liquid He temperatures and observed 
a diffuse maximum in the complex permittivity below 
40 K. The crystal structure of the pyrochlore PMN single 
crystal (Fd3m space group) was determined in detail 
by Wakiya et al£ High-frequency dielectric properties, 



including microwave, THz and infrared frequencies, have 
not been investigated, although they could be useful 
for understanding the reported diffused and frequency 
dependent maximum of the complex permittivity. The 
present report aims to fill this gap in the literature. 
We will show that our pyrochlore PMN ceramics does 
not undergo a diffuse phase transition (reported in 
Ref£), but a quantum paraelectric behavior for which 
the temperature dependent permittivity is caused only 
by anomalous polar phonons. In this way, it represents 
the first quantum paraelectrics with pyrochlore crystal 
structure. 



II. EXPERIMENTAL 

Pb1.83Mgo.29Nb!. 71 06.39 pyrochlore ceramic samples 
were produced by solid state reaction of mixed oxide 
powders described in details in Kefs/'- 7 . PbO (99.5%), 
Nb 2 5 (99.5%) and MgO (97%) powders were mixed 
and sintered at 880 °C for 8 hours. The pyrochlore cubic 
structure was verified by the X-ray diffraction. 

The dielectric response was investigated between 400 
Hz and 1 MHz from 10 K to 730 K using an impedance 
analyzer HP 4192A. The TEo n i composite dielectric res- 
onator method^ and network analyzer Agilent E8364B 
was used for microwave measurements at 8.8 GHz in 
100 - 350 K temperature interval. The cooling rate was 
2 K/min. 

Measurements at teraherz (THz) frequencies from 
7cm~ 1 to 33 cm" 1 (0.2 - 1.0 THz) were performed in the 
transmission mode using a time-domain THz spectrome- 
ter based on an amplified Ti - sapphire femtosecond laser 
system. Two ZnTe crystal plates were used to generate 
(by optic rectification) and to detect (by electro-optic 
sampling) the THz pulses. Both the transmitted field 
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amplitude and phase shift were simultaneously measured; 
this allows us to determine directly the complex dielectric 
response e*(w). An Optistat CF cryostat with thin my- 
lar windows (Oxford Inst.) was used for measurements 
down to 10 K. 

Infrared (IR) reflectivity spectra were obtained using 
a Fourier transform IR spectrometer Bruker IFS 113v 
in the frequency range of 20 - 3000 cm -1 (0.6 - 90 
THz) at room temperature, at lower temperatures the 
reduced spectral range up to 650 cm -1 only was stud- 
ied (transparency region of polyethylene windows in the 
cryostat). Pyroelectric deuterated triglicine sulfate de- 
tectors were used for the room temperature measure- 
ments, while more sensitive liquid-He-cooled (1.5 K) Si 
bolometer was used for the low-temperature measure- 
ments. Polished disk-shaped samples with a diameter 
of 8 mm and thickness of ^2 mm were investigated. 
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III. RESULTS AND DISCUSSIONS 

Temperature dependence of the real and imaginary 
parts of complex permittivity e* — e' — is" at various fre- 
quencies is plotted in Fig. [TJ One can see typical incipient 
ferroelectric behavior, i.e. increase in e' on cooling and its 
noticeble saturation at low temperatures. It is important 
to stress that within the accuracy of measurements no 
frequency dispersion of e' was observed between 400 Hz 
and 8.8 GHz at temperatures below 600 K. The small 
low-frequency dispersion above 600 K is caused by non- 
negligible conductivity of our sample. The pronounced 
e'(T) dependence is therefore caused by the softening of 
an excitation above 10 GHz. To reveal it we measured 
the THz dielectric spectra (see Fig. [2]) and IR reflectivity 
spectra (Fig. OH) below room temperature (RT). One can 
actually see very pronounced changes in the THz com- 
plex permittivity due to the polar phonon softening in 
the investigated range (see Fig. 03). 

In order to obtain all phonon parameters as a function 
of temperature, IR and THz spectra were fitted simulta- 
neously using the generalized-oscillator model with the 
factorized form of the complex permittivity^ 



J TOj ~ W + iuJ lTO] 



(1) 



where ojtOj an d ULOj denotes the transverse and longi- 
tudinal frequency of the j-th polar phonon, respectively, 
and 7tOj and jLOj denotes their corresponding damp- 
ing constants. e*(u) is related to the normal reflectivity 
R(w) by 



R(u) = 



(2) 



The high-frequency permittivity £oo resulting from 
the electron absorption processes was obtained from 
the room-temperature frequency-independent reflectivity 



FIG. 1: (color online) Temperature dependence of the real e' 
and imaginary e" part of complex permittivity in pyrochlore 
PMN ceramics at different frequencies. e(0) means the sum of 
phonon and electron contributions to the static permittivity, 
as obtained from the IR reflectivity and THz data fit. e" data 
are plotted only below 300 K, because its higher temperature 
values are influenced by the conductivity. Note the right scale 
fore". 



tails above the phonon frequencies and was assumed to 
be temperature independent. 

The real and imaginary parts of s*(lo) obtained from 
the fits to IR and THz spectra are shown in Fig. [4j Pa- 
rameters of the fits performed at 300 and 20 K are sum- 
marized in Table I. One can see a higher number of ob- 
served modes at 20 K than at RT. This is due to the 
reduced phonon damping at low temperatures, which al- 
lows to resolve a higher number of modes, which are prob- 
ably overlapping at RT. One can see in the e" (to) spectra 
of Fig. 02 (we remind that the frequencies of e" maxima 
roughly correspond to the phonon frequencies) that the 
most remarkable frequency shift with temperature is re- 
vealed by the lowest frequency mode below 30 cm , but 
also the mode near 130 cm -1 partially softens on cool- 
ing. The temperature dependence of the soft mode fre- 
quency wsm (left scale) and its dielectric strength Assm 
are shown in Fig. 03 Mainly this mode causes an increase 
in e' on cooling (see Fig. Q}- In the case of uncoupled 
phonons the oscillator strength fj — Aej.ui^Qj of each 
phonon is roughly temperature independent so that each 
softening of phonon frequency lotOj is connected with 
the increase of its dielectric strength Aej given by£ 
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In our case the soft-mode oscillator strength Jsm is tem- 
perature dependent, it increases twice from 3.3 .10 4 to 
6.5 .10 4 cm~ 2 . This indicates that the soft mode is cou- 
pled with some higher frequency mode assuming that 
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FIG. 2: (color online) THz dielectric spectra of the pyrochlore 
PMN at various temperatures. The soft mode frequency shifts 
down into the THz range on cooling, therefore the sample 
becomes less transparent at low temperatures (the noise in- 
creases) and the accessible spectral range narrows on cooling. 
Two frequency scales (THz and cm -1 ) are given. 
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FIG. 3: (color online) IR reflectivity spectra of the py- 
rochlore PMN ceramics at various temperatures below RT. 
Low-temperature spectra were obtained only below 650 cm" 
(transparency range of the polyethylene windows in the cryo- 
stat). 
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FIG. 4: (color online) Complex dielectric spectra from the 
fit to IR reflectivity and THz dielectric spectra at various 
temperatures. To see better the mode softening, e" (w) in inset 
the spectra are shown in the spectral range below 200 cm" 1 . 



fj — const. However, it is difficult to reveal with 
which mode is the soft mode coupled, because the oscil- 
lator strength of the most high-frequency modes is much 
higher than f$M and relatively small changes (below the 
limit of our fitting accuracy) in these modes could explain 
the increase of the fsM on cooling. 

The static permittivity obtained from the fit of the IR 
reflectivity is defined as 



e(0) 



(4) 



and is plotted in Fig. [T]in red solid dots. The values of 
e(0) are slightly lower (by about 20) than the experimen- 
tal values obtained at and below the microwave range. 
However, we believe that the disagreement is rather due 
to the experimental inaccuracy than due to real dielec- 
tric dispersion above the GHz range. So we conclude 
that the temperature dependence of the permittivity be- 
low 10 GHz in Fig.[T]is essentially due to anomalous polar 
phonons. 

The temperature dependence of the lowest-mode fre- 
quency below 35 cm" 1 was fitted with the Cochran law 



^SM 



(T) = y/A(T - T cr ) 



(5) 



where A is a constant and T cr is the critical softening 



4 



TABLE I: Parameters of the polar phonon modes in the pyrochlore PMN obtained from the fit of IR and THz spectra at 20 
and 300 K. Frequencies ujtOj, ^LOj an d dampings jroj, fLOj of modes are in cm , Aej is dimensionless, £00=5.87. 
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from the fit of the temperature dependence of permittiv- 
ity e(T), which (for classical paraclcctrics) should follow 
the Curie- Weiss law 

C 

± ± cw 

The result of the Curie- Weiss fit is shown in Fig. [6] 
in dashed line with the following fit parameters: 
e cwoo =43±1, Curie- Weiss constant C=(47800±500) K 
and critical temperature T cw =(-202±8) K. The Cochran 
critical temperature T cr =-285 K is somewhat lower than 
the Curie- Weiss critical temperature T cw , but if we con- 
sider that the extrapolated critical temperatures lie far 
below the investigated temperature range, the agreement 
between both values is reasonable. 

The Curie- Weiss fit in Fig. [5] deviates from the experi- 
mental data below ~ 50 K, because e'(T) levels off at low 
temperatures. Similar behavior was observed in incip- 
ient ferroelectrics as SrTi03, KTaC>3 or CaTi03 where 
the polar soft mode is also responsible for the observed 
e'(T)j 10 i 11 The soft mode does not soften completely, 
it levels off at low temperatures (typically below 30 K), 
mainly due to zero-temperature vibrations of light oxy- 
gen ions, which prevents the formation of long-range fer- 
roelectric order and permittivity divergence at low tem- 
peratures. Due to this phenomenon the incipient fer- 
roelectrics are also called quantum paraelectrics^ Note 
that pyrochlore PMN is the first quantum paraelectrics 
of pyrochlore crystal structure. 

The leveling-off of the low-temperature permittivity 
in incipient ferroelectrics was explained by Barretl^ al- 
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FIG. 5: (color online) Temperature dependence of the low- 
frequency soft mode usuQs& scale) and its dielectric strength 
(right scale). The dashed black line and solid red line show 
the results of the Cochran and Minaki fit of the soft-mode 
frequency, respectively (see the text). 



temperature. From the fit we obtained A = (1.86±0.01) 
cm" 2 FT 1 and T cr = (-285±11) K. So, pyrochlore PMN 
tends to the ferroelectric instability at negative tempera- 
tures. 

Theoretical critical temperature can be obtained also 
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FIG. 6: (color online) The temperature dependence of the 
experimental permittivity at 500 kHz and result of the Curie- 
Weiss fit (dashed green line) and the fit with Barrett formula 
(red solid line). Note the log temperature scale. 



ready in the beginning of 1950's. He derived the formula 



s'(T) = 



M 



^coth(|) 



IT i 



T 



£Bc 



(7) 



where M is constant, T% is the temperature below 
which the quantum fluctuations start to play a role 
(^fcsTi is the zero-point vibration energyM) and £boo 
marks the temperature independent part of the per- 
mittivity (this term was neglected in Ref.— , because it 
was very small in comparison to huge low-temperature 
e' in SrTi0 3 and KTa0 3 ). Our use of the the Bar- 
rett formula in Eq.( [7|) for fitting of e'(T) yields very 
good agreement with the experimental data (see Fig.[S]). 
We found M=(4.25±0.01)xl0 4 K, Ti = (96±8)K T =(- 
167±9) K and e Boo = 47±0.5. For T > T u \T X coth(g,) 
asymptotically approaches T and Eq.( [7]) becomes a 
Curie- Weiss law. Therefore also £boo — 47 is close to 
£cWoo =43. It is worth to note that the zero point vi- 
bration frequency \ksTi = lTHz=33cm _1 corresponds 
very well to the soft mode frequency (see Fig. [5|). Note 
also that the T\ parameter as well as the soft-mode fre- 
quency in the pyrochlore PMN qualitatively agree with 
analogous parameters obtained for SrTiOs and SrTi 18 03, 
although the values of permittivity in these materials are 
two orders of magnitude higher ) 13 ' 14 i 15 which is the con- 
sequence of different M parameters. 

In the only published dielectric data below 400 kHz 
Shrout and Swarta* observed similar essentially disper- 
sionless increase in e' on cooling down to 50 K as we 
did in Fig. [TJ We believe that the small dispersion be- 
low 50 K and small decrease in e' below ~ 30 K ob- 
served by Shrout and Swarta^ could be due to some 
defects (e.g. vacancies) in the crystal lattice of py- 
rochlore Pb1.s3Mgo.29Nb1.71O6. 39, which is substantially 
non-stoichiometric. 

It is clear that the soft mode frequency cannot follow 
the Cochran law (Eq. [S]) in the case of low-temperature 



quantum fluctuations. The correct low-temperature de- 
pendence of the soft mode frequency derived from the 
Barrett formula for permittivity can be found e.g. in the 
paper of Minaki et al. 16 



usm(T) 



cothf--) 
\2TJ 



T 



(8) 



where A is constant and Ti and To have the same 
meaning as in Eq.([7]). Note that Eq.([5} follows from 
Eq.([7]) and Lyddane-Sachs- Teller relation under assump- 
tion that the temperature dependence of static permit- 
tivity is caused just by softening of the soft TO mode. 
Result of the soft mode fit with Eq. ( |HJ) is shown by solid 
line in Fig. [5] where one can clearly see the leveling-off of 
the soft-mode frequency at low temperatures (unlike the 
Cochran fit). The fitting parameters are the following: 
A = (2.03±0.05) cm^K" 1 , Ti = (96±9)K and T = 
(-240±11)K. The fit parameters could be significantly 
improved if we would have more points in Fig. [5] espe- 
cially below 50 K. Nevertheless, one can see reasonable 
agreement of both Cochran and Barrett fits above 50 K 
as well as To and Ti parameters obtained from the Bar- 
rett fits of permittivity (Eq. [5]) and the Minaki model of 
the soft mode frequency (Eq. [5]). 

Let us compare the IR reflectivity spectra of pyrochlore 
PMN (Fig. E| with those of perovskite PMN. The lat- 
ter was first published by Burns and Dacol 3 together 
with the temperature dependence of the optical index 
of refraction n(T), which shows deviation from the lin- 
ear dependence below ~ 600 K. They explained the un- 
usual n (T) dependence by formation of polar nanore- 
gions. Surprisingly, the published IR spectrum of the per- 
ovskite sample 3 corresponds to our pyrochlore spectrum 
in Fig. [31 In later papers of other authors (Refsi 17 i 18 ' 19 i 20 ) 
the mutually similar (but different from Burns and Da- 
col's spectral) IR reflectivity spectra of perovskite PMN 
consisted of three distinct reflection bands typical for all 
cubic perovskite oxides. We stress that the infrared spec- 
tra in Refsj 17 ' 18 ! 19 ' 20 were obtained independently on ce- 
ramics, single crystals as well as on thin films. It ap- 
pears that the IR spectrum by Burns and Dacol 3 - be- 
longs to pyrochlore PMN. The rest of their data (n(T) 
and P(T)) were obviously obtained on perovskite PMN, 
because the peculiarities near the Burns temperature in 
the perovskite PMN were later confirmed in many exper- 
iments. 

It is of interest to compare the number of observed 
polar modes in the reflectivity spectra with the predic- 
tion of factor-group analysis: pyrochlore PMN crystal- 
lizes in Fd3m space group with 8 formula units per con- 
ventional unit cellpi i.e. 2 formula units per primitive 
unit cell. This means that on the whole, 66 lattice vibra- 
tional branches are expected. Pb ions are in 16d positions 
while Mg and Nb ions are in 16c positions^, both sites 
having D 3 d symmetry, while the O cations are in posi- 
tions 48f and 8b of C% v and T<j symmetry, respectively. 
The mode symmetries and their activities in IR and Ra- 
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man spectra can be obtained using standard tables^ with 
the following result for the T-point of the Brillouin zone 
(factor-group analysis) : 

IVdSm = 3A 2u (-) + 3E U {-) + 8F lu {x) + 4F 2u 
+4F 2g (xy, yz, xz) + A lg {x 2 , y 2 , z 2 ) 
+E g (x 2 + y 2 - 2z 2 , V3x 2 - VZy 2 ) + 2F lg . (9) 

It means that after subtraction of 1F± U acoustic mode, 
7 Fiu modes should be IR active, 4F 2g , lA lg and lE g 
should be Raman active, the rest of modes being silent. 
Table I shows that our fit of IR spectra required 17 
modes, much more than expected. The analysis in Eq. [S] 
assumes one effective ion in 16c positions instead of sta- 
tistically distributed Mg and Nb ions. Since the ions 
strongly differ in the mass, one could expect splitting of 
the modes in which these ions take part. If we take into 
account different Mg and Nb vibrations, the factor group 
analysis yields: 

T Fd3m = 4A 2 „(-) + AE U {-) + WF lu {x) + 5F 2u 
+4F 2g (xy,yz,xz) + A lg (x 2 , y 2 , z 2 ) 
+E g (x 2 + y 2 - 2z 2 , V3x 2 - W) + 2F lg . (10) 

In this case 9 polar F\ u modes are expected, which is still 
less than 17 modes observed at low temperature (see Ta- 
ble I). If we assume that the Pb cations and some of the 
oxygen anions are dynamically disordered among more 
equivalent positions with the avarage structure remain- 
ing cubic like in isostructural Bii.sZnNbi.sOy^, then 14 
modes could be IR active^ close to our experimental 
result. 

The problem with the excess IR active modes is also 
known from the perovskite PMN , where only AF\ U po- 
lar modes are allowed in Fm3m structure, although 7 
modes were observed^ The excess modes in the per- 
ovskite PMN were explained by polar clusters, which 
locally break the cubic symmetry into a rhombohedral 
one^ In the case of pyrochlore structure similar local 
symmetry breaking, if present, should probably be non- 
polar, because there is no indication for the existence 



of polar clusters. There is also no dielectric dispersion 
below the polar phonon range, in contrast to the PMN 
perovskite, where the huge dielectric dispersion appears 
just due to polar cluster dynamics J^- Therefore it ap- 
pears that new structural studies of pyrochlore PMN are 
needed to detect either a dynamical disorder of some 
atoms in the lattice or a non-cubic crystal symmetry. 



IV. CONCLUSION 

Our dielectric studies of pyrochlore PMN indicate 
quantum paraelectric behavior, i.e. the permittivity in- 
crease on cooling following the Barrett formula in the 
whole investigated temperature range. The permittivity 
shows no dispersion up to the microwave range and its 
temperature dependence can be explained by the soft- 
ening of polar optic modes. The zero-point vibrational 
energy ^UbTi =1 THz obtained from the Barrett for- 
mula corresponds very well to the soft mode frequency 
observed near 30 cm -1 . It is worth to note that the py- 
rochlore PMN is the first quantum paraelectrics with 
pyrochlore crystal structure. Our IR spectrum of py- 
rochlore PMN corresponds to IR spectrum by Burns and 
Dacol^, whose paper concerns the perovskite PMN. By 
comparing it with the later resultS ) 17 i 18 ' 19 i 20 it becomes 
clear that Burns and Dacol published (by mistake) the 
IR spectrum of the pyrochlore PMN. Since the IR experi- 
ment revealed more modes than expected from the factor- 
group analysis, we suggest that the structure should have 
at least locally lower symmetry than the cubic one. 
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